Quantum Criticality from in-situ Density Imaging 
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We perform large-scale Quantum Monte Carlo (QMC) simulations for strongly interacting bosons 
in a 2D optical lattice trap, and confirm an excellent agreement with the benchmarking in-situ den- 
sity measurements by the Chicago group [l| . We further present a general finite temperature phase 
diagram both for the uniform and the trapped systems, and demonstrate how the universal scaling 
properties near the superfiuid(SF)-to-Mott insulator(MI) transition can be observed by analysing the 
in-situ density profile. The characteristic temperature to find such quantum criticality is estimated 
to be of the order of the single-particle bandwidth, which should be achievable in the present or 
near future experiments. Finally, we examine the validity regime of the local fluctuation-dissipation 
theorem (FDT), which can be a used as a thermometry in the strongly interacting regime. 



In this recent decade, systems of ultracold atoms have 
become the most promising candidate for the realization 
of quantum simulators, which arc designed to explore 
various challenging and exotic many-body physics. One 
of the most well-understood ultracold systems is made of 
bosonic atoms loaded in an optical lattice, whereby the 
3D SF-MI phase transition had been observed in time-of- 
fhght(TOF) experiments [H and then directly compared 
to theory with great accuracy [|[. However, the TOF im- 
ages do not directly represent physical quantities of atoms 
inside the optical lattice and thus subject the quantita- 
tive study of critical phenomena to question, especially 
given the complexity of the expansion dynamics during 
the TOF imaging [jQ. As a result, several schemes 
have been proposed to determine the critical properties 
based on density-related quantities which have been 
successfully measured via in-situ spatially resolved den- 
sity imaging in a bosonic systems loaded into a 2D op- 
tical lattice [10. Before these theoretical schemes could 
be trusted and applied, an essential step is to validate 
them from the first principle calculation in the parame- 
ter regime of realistic experiments. 

In this paper, we perform an ab-initio QMC simula- 
tion [To| to have a direct comparison with the in-situ 
data of Chicago group [l| in a 2D optical lattice with 
a harmonic trap, and find an excellent agreement (see 
Figs, [lja) and (b)). We then calculate the finite temper- 
ature phase diagram for a uniform system, and justify 
the validity of the local density approximation (LDA), 
except near the phase boundary. We further quantita- 
tively confirm the universal scaling properties near the 
SF-MI quantum phase transition point through the fi- 
nite temperature density profile in a trapped system. 
The characteristic temperature, T*, to find such quan- 
tum criticality is comparable to the single particle band 
width, and therefore is achievable in the present or near 
future experiments. Finally, we show the regime where 
the local fluctuation-dissipation theorem (FDT) is valid, 
providing a useful thermometry in the strongly interact- 
ing regime. Our work therefore paves the way for the 
future experimental realization of quantum simulators. 
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FIG. 1: (Color online) Density profiles measured by Chicago 
group [l[ (blue squares) and by our QMC simulation (green 
lines) for (a) a shallow lattice and (b) a deep lattice. a(d at ) 
is the lattice constant (strength), and Er is the recoil energy, 
(c) is the finite temperature phase diagram of a uniform 2D 
system. Green thick line is the SF-MI phase boundary at T = 
0. Red(blue) lines are boundaries of SF(MI) at T/U = 0.01, 
0.05. 0.1, and 0.15 from left to right(from right to left). The 
regime of MI phase is defined in the text. The vertical dashed 
lines are the chemical potential range in Figs. [3] and [4] 



Our numerical simulation starts from the single-band 

Bose-Hubbard model: H = —JJ2(ij) '+ h.c.j + 

T Si ^*(^« — 1) — Ylti nijCijnii where hi and hi are the bo- 
son field operator and density operator at the ith lattice 
site. /u(r,) = // — F(rj) with /x the global chemical poten- 
tial and V(jCi) = Vt T x (xf + yf) being the isotropic trap- 
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FIG. 2: (Color online) (a) and (b) are typical density pro- 
files (red lines) and local compressibility (blue lines, arbitrary 
unit) in a trap system, obtained from the QMC calculation. 
Red dots and blue squares are results obtained by LDA, al- 
most identical to the QMC results. The thick green lines 
shows the local superfluid density within the LDA. Here N(V) 
stands for normal phase(vacuum). 



ping potential at lattice position, r.; = (xi,yi). Within 
a 5% uncertainty in the total particle number N, direct 
comparisons of the density images show excellent and 
quantitative agreements (sec Fig. QJa) and (b)) both in 
the weak and strong lattice strengths. A small deviation 
in the trap center of Fig. I[b) is due to some experi- 
mental equilibration issues [111 ]. Hence, our results con- 
vincingly indicate that the bosonic atoms loaded into an 
optical lattice is an excellent quantum simulator. 

In Fig. HJc), we show the finite temperature phase 
diagram in a uniform 2D system. The boundaries be- 
tween SF and normal phases (red lines) are determined 
by vanishing SF density[l2|]. Strictly speaking, MI phase 
doesnot exist at T > 0, but we can still definite its 
regime (blue lines) empirically by a small compressibil- 
ity, k = dn/dfi < 0.02/U. As a result, one can clearly 
see how both the SF and MI regimes shrink as T in- 
creases, leading to the normal phase in the intermedi- 
ate regime. Such finite temperature phase diagram of 
a uniform system can help in the determination of the 
phase boundaries in a trapped system within the LDA 
[f3| |. as shown in Figs. [Ha) and (b). The density profile 
(p{r)) as well as local compressibility (n(r) = dp(r)/dp) 
are obtained within the QMC simulation, but both of 
them can be well-reproduced within the LDA, except 
for the narrow regime near the phase boundary [14| . 
Such rounding of n(r) originate from the finite size ef- 
fect and the finite temperature Berezinskii-Kosterlitz- 
Thouless (BKT) transition. From our direct comparison, 
however, the phase boundary determined by the largest 
curvature points of k(t) is still very close to the position 
obtained within the LDA (within one or two lattice sites) , 
and also agrees with the onset of the local SF density. 
Therefore, practically speaking, one can still apply LDA 
to estimate the phase boundary even in the 2D trapped 
system at finite temperature. In the following, we will 
use another scheme to investigate the phase boundary as 
well as the quantum criticality near the SF-MI quantum 
phase transition from the density profile. 

It is well-known that near the second order phase tran- 
sition, the correlation length diverges and universal scal- 



ing behaviors are expected in many physical quantities. 
In the context of ultracold atoms, the universal scaling 
behaviors of a classical sup erfluid phase transition have 
been observed in a 3D |l5| | and 2D [l6| weakly interact- 
ing Bose gas when the temperature is near the transi- 
tion temperature. However, the scaling behavior near 
the quantum (zero temperature) phase transition point 
is still unexplored in current cold atom experiments. To 
provide useful indications for future experiments, we shall 
examine the finite-temperature calculation in quantita- 
tive detail, as proposed in Rcf. [6.7]. In the low-energy 
limit, it is well known that the SF-MI transition can be 
described by one of the following two effective models 
the dilute Bose gas model (DBG) for the phase 
transition in the edges of the MI lobe (for example, 
along the C\ — C*2 line of Fig. QJc)), and the 0(2) rotor 
model (or A"F-model) near the tip of MI lobe, where the 
particle-hole symmetry holds. 

Near the transition point, the free energy becomes sin- 
gular in an infinite system, and its singular part, f Sl 
should be a just function of p and T for a given J/U. 
As a result, a universal scaling theory gives / s (/i,T) ~ 
Ti+ 1 F((p-p c )/T 1 / zv ) [IG3, where d is dimensionality, 
z is the dynamical exponent, v is the correlation length 
exponent, /i c is the chemical potential at zero tempera- 
ture phase transition point, and F is the universal scaling 
function. As a result, we can obtain the singular part of 
the particle density, n s = n — n r = df s /dp,, and com- 
pressibility, = dn s /dp, as following (here 
n r (n r ) are the regular part of density (compressibility)): 



n s ((i,T) ~ Ti +1 -^G 



P~ Pc 



P - Mc 



T l/z 



(1) 



(2) 



where G(x) = dF(x)/dx and H(x) = dG(x)/dx are also 
universal functions. Combining them we could eliminate 
p — p c and obtain 



K s (n Sl T) = T* 



(3) 



where both n s and k s can be measured easily [T|. 

In the following analysis, we will first concentrate on 
the generic case, i.e the class of DBG model, and then 
show results of the 0(2) rotor model later. It is well- 
known that (z,v) = (2,1/2) and (z. v) = (1,1) respec- 
tively for them when fixing J/U (7llla|. According to Ref. 
Q, the cross point of n s /T at different temperature can 
give the value of \i c from Eq. (|T|), as clearly observed in 
our simulation data shown in Fig. [3Ja). Since we con- 
sider the regime with n r = 1 and n s > 0, it gives us the 
transition point C\ in Fig. [TJc). Similarly we can also 
obtain the position of C2 (not shown), so that their differ- 
ence gives the zero temperature gap for a given J/U. We 
also have performed the finite size scaling in a trapped 
system (proposed in Ref. [Iu3) to identify the precise 
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FIG. 3: (Color online) (a) n s /T as a function of fi(r) in 
a trapped system for various temperature, (b)-(d) are re- 
spectively the universal scaling function G, H, and H(G~ 1 ) 
as defines in Eqs. CD-©- Red(blue) lines are taken from 
the low(high) temperature data, 0.02 < T/U < 0.03(0.03 < 
T/U < 0.0625). The vertical and horizontal dotted lines in- 
dicate the critical values (Ci point of Fig. [He)). The rectan- 
gular window shows the parameter regime of (a). The sharp 
noise in (d) results from the fiat density at the cloud center. 
Here we use J/U = 0.03, V tr /U = 1.75 x 10" 4 , N ~ 2 x 10 4 , 
and exponents of the 2D DBG model here. 
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FIG. 4: (Color online) (a) and (b) are scaling function 
H(G^ 1 ) near the tip of MI lobe (J/U = 0.054), where the 
scaling exponents are from the DBG model and the 0(2) 
model respectively. Red(blue) lines are for the lower(higher) 
temperature as indicated in Fig. [3] (c) shows how k s changes 
as a function of temperature for several typical values of n s /T 
at J/U — 0.03 (same parameters as Fig. [3]in the DBG model). 
The vertical dashed line indicates the maximum temperature 
(T/U = 0.0625) used in Fig. [3] and the horizontal dotted line 
indicates k c at the critical point, (d) Similar to (c) but for 
J/U = 0.054 in the 0(2) rotor model. 



value of li c . The obtained result is identical (< 0.1%) to 
the result of a uniform system after finite size scaling [2(| > 
and also very close (< 1%) to the one we obtain above 
in a trap system (Fig. [3](a) ) . Our results varify the pro- 
posal to map out the zero temperature phase diagram in 
a realistic trapped system. 

After determination of /i c , in Fig. EJb) we show the 
rescaled singular density (n s /T) as a function of rescaled 
chemical potential, (li — li c )/T, in the DBG model, and 
find a very nice convergence for the universal scaling func- 
tion G (see Eq. ([T])) in a wider temperature regime. This 
is a clear evidence of the scaling theory near the critical 
point of SF-MI transition. Similarly, we also show the 
universal functions for H and H(G^ 1 ) in Figs. 02c) and 
(d) respectively. The later results are also shown for both 
the particle side and hole side of the phase transition 
points (i.e. the C\ and C*2 points in Fig. [TJ:)). 

In our calculation here, we find that the universal scal- 
ing function is very sensitive to the value of exponents, 
i.e. no intersection point for li c and no convergence of 
the scaling function can be obtained if using exponents 
of other universal classes. This can be also observed when 
we calculate the universal scaling function near the Mott 
tip, as shown in Fig [Ha) and (b). For these data at 
J/U = 0.054 (i.e. C 3 - C 4 line in Fig. Q^c)), exponents 
obtained from the 0(2) rotor model works much better 
than the DBG model, consistent with the general under- 
standing that the low energy effective model at the Mott 
tip is governed by the former. The crossover from the 
DBG model to the 0(2) rotor model near the Mott tip is 



of great interest, but is beyond the scope of this paper. 

From the experimental point of view, however, it is 
somehow more important to know the characteristic tem- 
perature, T* , below which the universal scaling behavior 
is observable. The effective model itself cannot provide 
this information, and therefore our ab-initio QMC data 
becomes very crucial. In Fig. Etc), we show how the 
universal scaling of the DBG model is affected by tem- 
perature at J/U = 0.03. We can see that k s keeps almost 
a constant in the low temperature limit up to some higher 
temperature, T* , above which their universal value starts 
to deviate due to higher order effects. For compressibil- 
ity near k c (i.e. /i ~ fi c ), such T* becomes much larger. 
If we set the minimum range for observing the scaling 
law in real space (within LDA) is about 10 lattice sites, 
the corresponding T* ~ 6 J is about the order of the sin- 
gle particle band width. Although such low temperature 
measurement is challenging, but should be still exper- 
imentally accessible in the near future. When J/U is 
close to the MI tip, due to the crossover from the DBG 
model to the 0(2) model, we find that T* is reduced and 
the range to observe the universal scaling behavior be- 
comes smaller (see Fig. [UJd)). Therefore, at least within 
the QMC simulation, our results confirm the possibility 
of observe the universal scaling theory near the critical 
point in the present experimental parameter regime. 

Accurate thermometry is the essence in the observa- 
tion of universal scaling behavior. Earlier scheme to fit 
the outskirt density profile with a noninteracting model 
hardly works in the strongly interacting regime we con- 
sider here. It has been proposed that the fluctuation- 
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FIG. 5: (Color online) (a) shows the calculated values of 
fcsTre(ri) (red dashed line) and (<5nf ) (green dotted line) and 
their ratio (x( r 0i blue solid line) for a typical trapped sys- 
tem, (b) and (c) show the contour plot of \ as a function of 
J/U and /j in a uniform 2D system at T/U — 0.05 and 0.25 
respectively. The solid and dashed lines in (a) indicate the 
SF and MI boundaries. The vertical line in (b) is the range 
of chemical potential in (a) . 



dissipation theorem (FDT) can be a feasible thermom- 
etr y s cheme (independent of any theoretical model) 
21 22j . but in realistic experiments, it is more practical 
to consider the local version of the FDT. We start from 
the identity, |^ 



{5N 2 )/k B T, where N = £\ n, is the 
total number operator. Therefore, the local compress- 
ibility can be related to the local density fluctuation by: 
K{vi) = = {5n})/k B T + Y, j ^fi i 8n j )/k B T, where 
the last term is the non-local density correlations. In 
the normal phase regime away from the critical point, 
this term should be negligible due to the short correla- 
tion length, i.e. k B Tn(Yi) = (8nf). In order to investi- 
gate the regime when the local FDT is valid, we define 
x(r,) = k B T 'ft(rj) / (Snf) , and show the result of a typical 
parameter in Fig. Eta), where there is a wide regime (not 
only in the outskirt regime) to find a constant plateau, 
which can estimate the system temperature with a rea- 



sonably high precision. More generally, in Fig. [5jb) and 
(c), we find the regime for % « 1 is alway in the normal 
phase regime of small J/U (and away from the Mott tip), 
and then occupies the entire lower-left corner when the 
temperature is higher. Our results show that even is a 
pretty low temperature (but still T > J as in the present 
experimental limit), the wide range of normal phase helps 
the justification of the local FDT as a thermometry. This 
method takes the advantage that only local data mea- 
surement is needed, and can be also used to study the 
non-local correlation effects in the current experiments. 
When deep inside the SF regime (more particle number 
per site) or when T < J, non-local correlations become 
significant and the local FDT should be corrected with 
higher order terms as suggested in Ref. [2lj |. 

In conclusion, we show the direct comparison between 
the in- situ experimental data and the ab-initio numerical 
calculation for bosonic atoms inside a two-dimensional 
optical lattice with a trapping potential. We map out the 
quantum phase diagrams in a uniform system by using 
the knowledge of the density profile only. Besides, we 
note that the universal scaling function as well as other 
quantum critical properties can be observed if the system 
temperature is lower than the single particle band width. 
Finally we investigate the regime when the LDA and the 
local FDT works, where one can use them to extract 
phase diagram and system temperature. Our results can 
be directly applied to the current experimental and for 
future investigation of the quantum simulator. 
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